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The Jaynes-Cummings model without the rotating-wave approximation can be solved exactly 
by extended Swain's ansatz with the conserved parity. The analytical approximations are then 
performed at different levels. The well-known rotating-wave approximation is naturally covered in 
the present zero and first approximations. The effect of the counter rotating-wave term emerges 
clearly in the second order approximation. The concise analytical expressions are given explicitly 
and can be applicable up to the ultra-strong coupling regime. The preliminary application to the 
vacuum Rabi splitting is shown to be very successful. 
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I. INTRODUCTION 

The interaction of light with matter is a fundamen- 
tal one in optical physics. The simplest paradigm is a 
two-level atom coupled to the electromagnetic mode of 
a cavity. In the strong coupling regime where the cou- 
pling strength g/ui (a; is the cavity frequency) between 
the atom and the cavity mode exceeds the loss rates, 
the atom and the cavity can repeatedly exchange excita- 
tions before coherence is lost. The Rabi oscillations can 
be observed in this strong coupling atom-cavity system, 
which is usually called as cavity quantum electrodynam- 
ics (QED) [ij. Typically, the coupling strength in cavity 
QED reaches g/ui ~ 10^^. It can be described by the 
well-known Jaynes-Cummings (JC) modelQ. 

Recently, for superconducting qubits, a one- 
dimensional (ID) transmission line resonator [31 or 
a LC circuit [J-|6| can play a role of the cavity, which 
is known today as circuit QED. More recently, LC 
resonator inductively coupled to a superconducting 
qubit|8l-[l0l| has been realized experimentally. The 
qubit-resonator coupling has been strengthened from 
g/ui ~ 10"'^ in the earlier realization^, a few percentage 
later [a, 0|j to most recent ten percentages |8l-[l0{. In 
cavity QED system, the rotating-wave approximation 
(RWA) is usually made, however, in the circuit QED, 
due to the ultra-strong coupling strength g/w ~ 0.1 
[8l-[l0| , evidence for the breakdown of the RWA has been 
providedfs]. Therefore, counter rotating-wave terms 
(CRTs) in the JC model should be included. Recently, 
many works have been devoted to this qubit-oscillator 
system in the ultra-strong coupling regimelll|? 4la|. 

Actually, the JC model without the RWA has been 
studied extensively for more than 40 years. A incomplete 
list is given by Refs. |16l - l29| . Two main schemes are 
employed. One is based on the photonic Fock states [Ig- 
[2]| with the pioneer work by Swain[la|. The continued 
fractions are present and the solution ten becomes very 



intricate. The other is on the basis of various polaron-like 
transformations or shifted operators, which are basically 
photonic coherent states approaches |22h29| . The very 
accurate solution can be obtained readily, but the infinite 
photonic Fock states should be involved. 

The RWA eigenstates only include two bare states, 
which have facilitated earlier clean investigations on var- 
ious quantum phenomena in quantum optics. Surpris- 
ingly, one or a few dominant terms in the eigenstates of 
the JC Hamiltonian beyond the RWA ones are still lack- 
ing or not given explicitly until now, to the best of our 
knowledge. What emerges first beyond the RWA results 
may be very useful to analyze the effect of CRTs on var- 
ious phenomena at the microscopic level. In this sense, 
a few dominant terms are more helpful than the exact 
solution including infinite bare states. 

In this paper, by using the conserved parity, we extend 
Swain's wavefunction to the JC model without the RWA. 
We will not follow the usual exact diagonalization rou- 
tine. Alternatively, we derive a polynomial equation with 
only single variable, just the eigenvalue that we seek. The 
solutions to this polynomial equation can give exactly all 
eigenfunctions and eigenvalues for arbitrary parameters. 
Moveover, we can perform approximations step by step 
with the help of these exact solutions. The zero and first 
order approximations will exactly recover the RWA re- 
sults. The dominant effect of the CRTs emerges in the 
second order approximation. 

Without the RWA, the Hamiltonian of a two-level atom 
(qubit) with transition frequency A interacting with a 
single-mode quantized cavity of frequency a; is 

H ^ —(T^+uja^a + g{a^ +a)crx, (1) 

where g is coupling strength. Ox and ciz are Pauli spin- 
1/2 operators, a' and a are the creation and annihilation 
operators for the quantized field. Here, (5 = A — a; is 
defined as the dimensionless detuning parameter. The 



energy scale is set w = 1 here. 

The RWA is made by neglecting the CRTs, a^a+ + 
acr_, then one can easily diagonalize the Hamiltonian 
and obtain the eigenfunctions as[30j 
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For later use, we also list the relevant eigenvalues 
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where i?„ — ^(5^ + 4^^ (n + 1). In the ground state 
(GS), the qubit is in GS and the photon is in a vacuum 
state. The GS energy is ^° = -f . 

Associated with JC Hamiltonian with and without the 
RWA is a conserved parity 11 , such that [H, 11] = 0, which 
is given by 
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CTajCXp 



(iTlN\ , 



(5) 



where N = a^a is the bosonic number operator. 11 has 
two eigenvalues ±1, depending on whether the excitation 
number is even or odd. The above two states (6) and (7) 
with even n are of odd parity and with odd n even parity. 
The ground state is of even parity. The RWA results for 



the first 8 energy levels at resonance, i?i(2)n 



g\/n + 1), are given in Fig. [T] (a) for later comparisons. 



II. EXACT SOLUTION WITHOUT THE RWA 

First we introduce a scheme to obtain the exact solu- 
tions to the JC model without the RWA. For convenience, 
we can write a transformed Hamiltonian with a rotation 
around an y axis by an angle ^. In the Matrix form it 
takes 



H 



a^a + g{a^ + a) 
_A 

2 



A 
' 2 



''a — g{a^ + a) 



(6) 



About 40 years ago, Swain first proposed an ansatz for 
the wavefunction in the photonic Fock states[l6J, which is 
also a starting point of the standard numerically exact di- 
agonalization (ED) scheme. Since then, various methods 
have been developed along this line [17l - l20l | , but the con- 
served parity has not been considered, to our knowledge. 
Therefore continued fractions are unavoidably presented 
in the expressions for the eigensolutions. 

We also proceed along this line, but the parity is in- 
corporated in the Swain's ansatz, which is given by 
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with -l-(— ) stands for even (odd) parity, M is the trun- 
cated number. The Schrbdinger equation gives 

2J a^ac„ |n) 4- c/^J "^" (\/n |n - 1) + Vn+ l\n+l)) 



A 



T^E(-1)"^"I") = ^^"l")- 



(8) 
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Left multiplying the photonic states (m| gives 



mcm + gVm + lc™+i + g^/mcrn-i T y (-l)™Cm = ECm, 
then we have a recurrence relation 
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By careful inspection of Eq. ([T]) , one can find that cq is 
flexible in the Schrbdinger equation where the normaliza- 
tion for the eigenfunction is not necessary, so we select 
Co = 1.0. Then we have 
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Once the first two terms are fixed, the coefficients of the 
other terms higher than a+ should be determined by the 
recurrence relation Eq. ^ 

For 771 — M, the terms higher then (a"*") are ne- 



± glected, we may set cm+i = 0, then we have 



2 ^ ' 



CM - gVMcM-i = 0. (10) 



Note that this is actually a one- variable polynomial equa- 
tion of degree M. The variable is just the eigenvalue we 
want to obtain. It is expected that the roots of Eq. (fTO|) 
would give the exact solutions to the JC model without 
the RWA if M is large enough. 

To obtain the true exact results, in principle, the trun- 
cated number M should be taken to infinity. Fortunately, 
it is not necessary. It is found that finite terms in state 
d?]) are sufficient to give exact results in the whole cou- 
pling range. The usual numerical exact diagonalization 
can readily give the energy levels and their wavefunctions 
in the JC model, which can be regarded as a benchmark. 
We will increase the truncated number M until the rel- 
ative difference of the energies obtained from the roots 
and the standard value is just less than 10"^, which sets 
the criterion for the convergence achieved. Interestingly, 
for coupling constant 5 < 0.1 for three typical atom 
frequency A = 0.1,1, and 1.5, the truncated number 
M = 15 in the polynomial equation can give exactly 20 
energy levels by the above criterion for convergence. For 
g = 1.0 and 2.0, the polynomial equations with M = 25 
and 40 can yield about 20 energy levels exactly. In fact, 
all above calculations can be immediately done in an or- 
dinary PC. 
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Transforming back to the original frame gives 

a|ff,0), 






co|0) 
co|0) 



Interestingly, the first element in the zero-order approx- 
imation actually gives exactly the ground-state in the 
RWA. The other solutions are g irrelevant and therefore 
omitted. 



FIG. 1: The first 8 energy levels as a function of coupling 
constant g at resonance. The exact results without the RWA 
are presented by black solid line, (a) The RWA results (red 
dashed line) and (b) The analytical results in Eqs. (|26p . (|27|l . 
(I30|l , and (|32|l of the second-order approximation (red dashed 
line) . 



For late use, the first 8 exact energy levels as a function 
of the coupling constant g obtained by the above poly- 
nomial equations with M = 25 at resonate case (A = 1) 
are presented in Fig. 1(a) and (b) by black solid lines. 
The parity is not changed after the level crossing. 

We then try to follow the energy curves to get some 
analytical approximate results in the next section. 



III. ANALYTICAL RESULTS WITHOUT RWA 

The recurrence relation Eq. (jH]) can be simplified to a 
tridiagonal form 



fno{E) 1 
1 ni{E) V2 

\/5 1^2 (-B) 



\ 



\0 

where 







^ Om(^) / 



/CO ^^ 

Cl 
C2 



0, 



^m(-E) — - 

5 



m-£;Ty(-ir 



(11) 



The eigensolutions can be obtained from the zeros of the 
determinant of the Matrix, 

The zero-order approximation f2„i (E) = gives 

£;" = mT(-l)" — • 



The lowest energy is 
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;-0 



The ground-state is of even parity, so the state is 
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(12) 



A. The first approximation 

The first approximation is made by selecting the ma- 
trices with two order along the diagonal line 



\/m+T f7„,(S) 






(13) 



It is expected that the m ^ th two order determinant 
would contain the information of the two levels with same 
parity. Comparing with the RWA results, it can be in- 
ferred that even m is corresponding to odd parity and 
odd m to the even parity. Fortunately, we have the 
same n{E) for any value of m 



n^{E) = - 
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Then we have following quadratic equation 

A\ / . „ A' 
"2 



m- E+ — ] { m + 1 - E - — ] - {m + 1) g^ = 0, 



which yields the eigenvalues 
E^ri 



l±yil-Af+4g^m + l). 



(14) 



They are just the RWA ones in Eqs. ^ and (g]) . 

According to the wavefunction Eq. ([7]), the eigenstate 
is then obtained as 



H' = (^ 



Cm \m) +c„i+i \m + 1) 
±(-1)™ [cm \m) - Cm+l \m 
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(15) 



By the above relation between m and parity, we always 
have ±(—1)™ = —1. Transforming back to the original 
frame gives 
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which are just the eigenstates under the RWA in Eq. ([2]) 
for excited states. 

So in the first approximation, we can not obtained re- 
sults more than the RWA ones for all excited states. The 
effect of the CRTS should only emerges beyond the first 
approximation. 



B. The second order approximation 



with 



Naturally, the second order approximation is per- 
formed by reducing to the m—th three order determinant 



y/m+l n ,n+i{E ) y/m + 2 

v^^rT2 n„^+2{E) 



0. (17) 



For more concise, we only consider the resonant case S ~ 
0. It is straightforward to extend to finite detunings. 

Two univariate cubic equations for even and odd par- 
ity can be explicitly derived for any m — th three order 
determinant. Three roots yi < 2/2 < 2/3 for each univari- 
ate cubic equation can be obtained easily by the formula 
presented in the Appendix A. Comparing with the exact 
ones obtained above, one can find that the second roots 
2/2 are the true solutions. Especially, the first root 2/1 for 
m = with even (odd) parity gives the energy in the 
GS (the first excited state). With these results in mind, 
the general solutions can be grouped as the following two 
cases, and all eigenvalues and eigenfunctions can be given 
analytically in the unified way. 



1. m — 2k with even parity and m = 2fe + 1 wtth odd parity 

For both m = 2k with even parity and m = 2k + l with 
odd parity (fc = 0, 1, 2, ... ), we have the same univariate 
cubic equation in the following form 
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3m+-] E"^ 



8 - (9m + 27) .g2 
^J[4+{6m + 9)g^ 

Especially, m — with even parity will give the GS ad- 
ditionally. The GS energy is given by the first root 2/1 
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I + v/4 + 9.g2 [cos + V3 sin 9] 
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with 



= - arccos 
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The states in this case all takes the form 
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(21) 



2. m = 2k + 1 with even parity and m — 2k with odd parity 

For both m — 2A:+1 with even parity and m — 2k with 
odd parity, we have the same univariate cubic equation 
in the other form 
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According to the Appendix A, we have 

A = 4+(6m + 9).g^ 

B = -8m - 12 - {l2m'^ + 22m - 3) g'^ 



C = (2m + 3)-' / + [6m^ + 13m^ + |"^ - f ) 3^ 



^2 „4 



-|-4m -I- 12m -I- 9. 

It can be readily proven T — B^ — 4AC < in this case, 
so there are three different real roots. Note above that 
the energy level is given by the second root 2/2, so 



+ 2m2 + 6m, + 



0. 

(22) 



Similarly, we have 

A = (6m + 9) 5^ + 4, 

B = (-12m2-38m-2l).g2 + (_8m-4), 

C = (-2m-3fg^+i6m^ + 29m^ + —m+ — 
\ 2 4 

-|-4m^ + 4m + 1. 

One can also readily prove T = B^ — 4 AC < 0, So there 
are also three different real roots. The energy level is 
given by the second root 2/2 , 



E = 



with 



(3m +1) + v/(6m + 9).g2 + 4 [cos 6 - VSsi 



(23) 
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(3m + I) + ^4 + (6m + 9) g^ [cos - \/3 sin O] 
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Especially, m = with odd parity will yield the first 
excited state additionally. The corresponding eigenen- 
ergy is given by the first root yi 



El = yi 
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with 



= - arccos 
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The states in this case all takes the form 
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3. Unified expressions 

For the future use, we will give the unified expression 
of the eigenvalues and eigenfucntions, which are corre- 
sponding to those in the RWA one by one in the follow- 
ing. 

Set TO = 71 — 1 in II B(l) and to = n in II B(2), the 
eigenvalues in Eqs. (|19p and ((23)) can be expanded in 
terms of g as 

1 
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The corresponding eigenstates take the form 

l"'>«((c',M|»-l> + c.«|„+l))). (^«) 

c„ \n) + Cn+2 \n + 2) 
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It should be pointed out that the GS state and the first 
excited state can not be brought into the above general 
expression for n = 0. The GS energy and the GS state 
with even parity are 
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and the first excited state with odd parity are 
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In any case, the ratios of coefficients in the second ap- 
proximation are 



■ Cj7i-|_i . Cn 
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m+2 
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where the sign =F in Eq. ([Tl]) for ftmiE) for any eigen- 
states is only parity dependent. 

It is interesting to note that the first two terms in Eqs. 
dMl), (EH) are no other than the RWA results by Eqs. §\) 
and dl]) at resonance. The additional terms appear just 
because of the presence of the CRTs. Also the eigenfunc- 
tions in Eqs. ((28|) and (|29p contain the components of the 
RWA ones in Eq. ([2]). The other bare state which can not 
be generated by the rotating-wave terms also emerges. 
This is just our answer to the question presented in the 
title of this paper. 

The analytical results for the energy levels in the sec- 
ond order approximation are collected in Fig. [IJb) with 
red lines. It is shown that for g < 0.4, the present second 
order approximation can give reasonable good results. 
More over, it should be deeply impressed that analytical 
expressions are almost exact for remarkable wide cou- 
pling regime 5 < 0.2. So it could become a solid and 
concise platform to discuss the effect of CRTs on various 
physical phenomena in the present experimentally acces- 
sible systems. Note that the present maximum value for 
the coupling strength in the superconducting qubit cou- 
pled to a circuit resonant Q has reached g — 0.12, to our 
knowledge. The application to the Vacuum Rabi split- 
ting is performed in the next section as a first example. 



IV. VACUUM RABI SPLICINGS 

If we pump the dressed atom from its ground to an 
excited state, it will decay to the dressed ground state 
through spontaneous emissions. Under the RWA, when 
the atom is excited by the operator, V — \e){g\ + \g){e\, 
from the ground state \g,0 > , the emission spectrum 
has two peaks with equal height. The distance of the 
two peaks, 2g, is just the vacuum Rabi splitting ^I, 3^. 

Without the RWA, we have two choices for the initial 
states. When the CRTs are included, the photon in the 
ground state are no longer a vacuum state, as shown in 
Eq. pip . In the framework of the second order approxi- 
mation, we first use V to excite the atom from the ground 
state at resonance 



\VRS)'^^^ ^V\GS) 
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which can be expanded in terms of the eigenstates with 
odd parity. Note that only the following 4 excited states 
are included 

\VRS)^^^ = vi \1EX)+V2 \2EX)+V5 \5EX)+ve \6EX) . 

(35) 
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We can derive the two peaks up to o(g^) 
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FIG. 2: (Color online) The peak heights (dashed lines) in 
emission spectrum as a function of coupling constant at res- 
onance for the initial states (a) V\GS) and (b)|e,0 >. The 
exact results are marked by black line, the present analytical 
results by dashed lines. 



The probabilities hi — {h'il can be expressed in g as 
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The atom will decay from the initial state to the dressed 
ground state with an emission spectrum. The heights 
of the peaks in the spectrum are proportional to the 
square of the probability of the corresponding eigen- 
states. Therefore one may find two main peaks, and the 
other two peaks are too small to be visible, by the above 
4 probabilities. 

In Fig. [5] (a), we plot the peak heights from the first 
two excited states hi and /i2 with dashed lines. The ex- 
act numerical results for the heights are also list with 
the solid lines. It is interesting to note that the present 
analytical results for the main peak height agree excel- 
lently with the exact ones in a wide coupling regime 
(0<5<0.2). 

With the increase of the coupling strength, the third 
peak from the sixth excited state jhA becomes visible. 
Recent full numerically exact study [33] showed three 
peaks (not four peaks) ai g = 0.8 in their Fig. [IJb). 
Our analytical results are consistent with this exact study 
qualitatively. If the third peak is visible, the coupling 
constant should exceed 0.5, the present second order ap- 
proximation can only describe it qualitatively. 

The other initial states is usual one |e, > which only 
include the first and second excited states in the frame- 
work of the second-order approximation 
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Those peak heights are also list in Fig. [2jb) with red 
lines. It is also shown that the analytical results in this 
case is consistent perfectly with the exact ones in a wide 
coupling regime (0 < g < 0.2). 

In both initial states, the level difference for the first 
two excited states will give the vacuum Rabi splitting 
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which is smaller than the RWA one by a small amount 

For the recent experimentally accessible ultra-strong 
coupling constant g = 0.1, the effect of CRTs on the 
vacuum Rabi splitting only results in a tiny shift around 
0.0025(aj), which is too small to be distinguished in the 
experimental data. While the ratios of the two heights 
hi/h2 for the first and the second initial states can be 
evaluated as 1.107 and 0.905 respectively, which are how- 
ever large enough to be seen experimentally. 

Finally, we would like to give some remarks. As shown 
in Fig. [U in a wide coupling regime {0 < g < 0.2), 
the difference between the RWA energy and the present 
second-order approximate energy is very subtle, but the 
states are essentially different. Some bare states in the 
latter are absent in the former. This is also the reason 
that the difference in vacuum Rabi splitting is invisible, 
but in the peak heights is evident. In the JC system, 
the accuracy of the eigenenergy is easy to ensure within 
various approaches, but it is not so crucial, in our opinion. 
The components in the eigenstates are very important, 
and play the dominate role in many physics processes. 



V. SUMMARY 

In this paper, the JC model without the RWA is 
mapped to a polynomial equation with a single variable, 
the eigenvalue, by the bosonic Fock space and parity sym- 
metry. The solutions to this polynomial equation recover 
exactly all eigenvalues and eigenfunctions of the model 
for all coupling strengths and detunings. Furthermore, 
the analytical results are presented at different stages. 
The first approximation in the present formalism repro- 
duces exactly the RWA results. The effect of the CRT 
emerges clearly just in the second order approximation. 
All eigenvalues and eigenfunctions are derived analyti- 
cally. It is shown that they play dominant role up to 
the remarkable coupling strength g — 0.2, suggesting 
that they could be convincingly applied to recent cir- 
cuit quantum electrodynamic systems operating in the 
ultra-strong coupling regime up to g = 0.12. Applica- 
tions of analytical results to the vacuum Rabi splitting 



are performed. Different lieiglits of the two main peaks 
are given explicitly, which agree well with the exact ones 
in a wide coupling regime. The concise analytical results 
only including three bare states will be very useful for 
the exploration of effects of CRTs on various phenomena 
in the ultra-strong coupling regime. 
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Appendix A: Solutions to univariate cubic equation 

The univariate cubic equation can be generally ex- 
pressed as 

x^ + bx^ + C.T + rf = 0. 



Its solutions can be found in any Mathematics manual. 
If 

r = B^ - AAC < 0, 

with 

A = b^ - 3c, B = be - 9d,C ^ c^ - 3bd, 

there are three different real roots with yi < j/2 < 2/3 
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